Localized Spins on Graphene 
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The problem of a magnetic impurity, atomic or molecular, absorbed on top of a carbon atom 
in otherwise clean graphene is studied using the numerical renormalization group. The spectral, 
thermodynamic, and scattering properties of the impurity are described in detail. In the presence 
of a small magnetic field, the low energy electronic features of graphene make possible to inject spin 
polarized currents through the impurity using a scanning tunneling microscope (STM) . Furthermore, 
the impurity scattering becomes strongly spin dependent and for a finite impurity concentration it 
leads to spin polarized bulk currents and a large magnetoresistance. In gated graphene the impurity 
spin is Kondo screened at low temperatures. However, at temperatures larger than the Kondo 
temperature, the anomalous magnetotransport properties are recovered. 

PACS numbers: 73.20.Hb, 73.20.-r, 73.23.-b, 81.05.Uw, 99.10.Fg 



Graphene is a two dimensional material made of car- 
bon atoms arranged in a hexagonal lattice. Its structural 
stability and unusual electronic properties P, S S 0) S [1] 
make it an excellent candidate for technological applica- 
tions. The low energy electronic structure corresponds 
to massless, chiral, fermionic quasiparticles described by 
the Dirac equation. Graphene is a semimetal that can 
be globally or locally doped with electrons or holes using 
gate electrodes. These characteristics triggered an in- 
tense activity that ranges from the search of new devices 
to the study of new scenarios for Dirac fermions 0, Q . 
Important advances have been made in the preparation 
and characterization of this material. One of the ongo- 
ing goals is to incorporate spintronic effects in graphene 
and this requires the development of simple tools for the 
manipulation and control of the carrier's spins. There 
are already some advances in this direction like the in- 
jection of a spin polarized current from ferromagnetic 
electrodes 0, S]. There are also some theoretical pro- 
posals involving the use of edge states to transport spin 
polarized currents Despite the intense activity 

in the area, the properties of graphene with magnetic 
impurities, atoms or molecules, have received less atten- 
tion. Previous works using mean field approaches already 
pointed out the unusual behavior of some properties as 
well as the possibility of controlling the magnetic struc- 
ture of the impurity with electric fields EE El. The 
magnetic screening of an impurity spin, the Kondo effect, 
in systems with graphenelike pseudogaps has also been 



analyzed by several authors IJ, 3, 1 



In this work we address the problem of graphene with 
magnetic impurities and show that the peculiar electronic 
properties of this system lead to some interesting new ef- 
fects. In particular, we show the potential use of these 
impurities to inject and generate spin polarized currents. 
When an impurity is adsorbed on top of a carbon atom, 
the impurity levels acquire a finite lifetime, that is, the 
spectral function shows broad peaks. In the undoped 



case, despite of the broadening of its levels the impu- 
rity behaves as a free spin at low temperatures and we 
show that with a small magnetic field, such that the Zee- 
man energy is larger than the thermal energy fc^T, a 
non-magnetic STM tip can be used to inject spin polar- 
ized electrons with an extraordinary efficiency. In this 
regime we also show that the bulk transport properties 
present interesting features: in the absence of electron- 
hole symmetry and with a finite impurity concentration 
there is a large magnetoresistance and the transport cur- 
rent is spin polarized. Conversely, a magnetic impurity 
in doped graphene leads to the Kondo effect. While in 
general the Kondo temperature Tk is small due to the 
small density of states at the Fermi energy of graphene, 
in some cases it could be well above the experimentally 
accessible limits. 

Our starting point is the Anderson model describ- 
ing an impurity with a single orbital of energy Sd and 
Coulomb repulsion U hybridized with the conduction 
electron states with a matrix element Vhyb. The Hamil- 
tonian of the system is then H = iJimp + -f^graph + -ffhyb, 
where the first term is given by 



(erf - a^iB)dld^ 



(1) 



Here dj^ creates an electron with spin a at the impurity 
state and /iB is the Zeeman energy shift due to an exter- 
nal in-plane magnetic field B. The Hamiltonian of the 
graphene layer is 

iJgraph = -tJ^^if^Wk^bk. + r{k)bl^a^,, (2) 



k.a 



where flj^^ and 6j^^ create electrons with spin a and 
wavevector k on sublattices A and B, respectively. The 
hopping matrix element t is of the order of 2.7eV @ 
and (j){k) — ^jC^^'^^ with {Sj} the three vectors con- 
necting one site with its nearest neighbors. As a re- 
sult there are two bands of width 3t that touch each 
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and the energy difference between the peaks is smaller 
than U . This is due to the interpay between a Hartree 
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FIG. 1: (Color online) Color map of the spectral density of 
a magnetic impurity on graphene. Parameters are U — 2eV, 
^hyb = leV. (b) STM spectra at the impurity for tc/td ~ 
(solid lines) and tc/td = 0.3 (dashed lines). Here Sd ~ — leV. 
(c) Same as (b) with Ed = — 0.5eV and tc/td ~ 0, 1.0. 



other at the corners of the Brillouin zone (Dirac point). 
Finally, assuming that the impurity is adsorbed on a 
site of sublattice A, the hybridization Hamiltonian is 
Hhyh = ^^J2ka^ia'^<y + ^^l^^ka , with N the number 
of unit cells in the sample. 



One of the main ingredients that determines the nature 
of the solution is the local density of states of the host 
material. Close to the Dirac point (i? = 0), the local den- 
sity of states, per atom and spin, can be approximated by 
p{E) = a \E\ with a — Auc/inti^Vp, where Auc is the unit 
cell area and is the Fermi velocity. In what follows we 
use this form with a high energy cutoff D. The other rel- 
evant parameters (e^, U, Vhyh) depend on the nature of 
the impurity, in particular U is of the order of a few eVs 
for transition metal impurities and smaller for molecules. 
We solve the problem using the extensions of Wilson's 
numerical renormalization group (NRG) method that al- 
low to describe a non constant density of states for the 
host material and to improve the accuracy of the high 
energy features 15, l3|- We have studied this model nu- 
merically for a wide range of parameters. In what follows 
we will focus on the localized spin regime where the av- 
erage number of electrons in the impurity level is of the 
order of one. 

In Fig. [TJa) we present a color map of the impurity 
spectral density A{u}) for the undoped case- the Fermi 
energy Ep laying at the Dirac point. The maximums of 
A{uj) are shifted from the bare parameters £d and Sd + U 
shown in the figure with dashed lines. While for a general 
p{E) some shifts are expected, here the shifts are large 



correction and the hybridization self-energy Figs, 
lb and Ic show the impurity spectral functions for two 
different values of Sd- Note the absence of a Kondo peak 
at Ep. These spectral densities could be measured with 
a STM where electrons from the microscope tip tunnel 
to the impurity sensing the local density of states. If the 
resonant level is close to the Dirac point {\ed\ < leV) the 
STM can measure the resonance. In general there is some 
leaking of electrons that tunnel to the substrate generat- 
ing Fano structures. The STM differential conductance 
at low temperatures is then given by [l8| 



G(14) 



4e^ 
Trfi 



rptAsTAiiEp + eVb). 



(3) 



where V^, is the voltage drop from the tip to the sample, 
and pt is the tip density of states at the Fermi energy. 
The quantity Astm (E) is the spectral function of the op- 
erator (icV'i + tddlr)/i, tc and td are matrix elements for 
the tunneling of an electron from the tip to the conduc- 
tion band states and to the impurity orbital, respectively, 
t = it1 + f^y/^ and is the field operator that creates 
an electron in a graphene state centered below the tip. 
In what follows we consider that the tip is on top of the 
impurity and for simplicity we take i/'J. = N~i Sfco- '^lo-- 
Even for tc — td the effect of tc is very small due to the 
smallness of p{E ^ Ep), see Fig. [1] and the STM gives 
direct information of the impurity spectral density. 

In the presence of an external magnetic field the impu- 
rity is polarized and the spectral densities become spin 
dependent. In Fig. [2] we present results obtained at low 
temperatures and low fields. Almost all the weight of the 
spin-resolved spectral densities is transferred to a single 
peak at the renormalized energies (for the spin up) and 
£d + U (for the spin down). For the small magnetic field 
used in the calculation the Zeeman shifts of the peaks 
are not appreciable. These results suggest that the spin 
dependent STM differential conductance at high voltages 
becomes very different for the two spin orientations. To 
estimate the current of spin-cr electrons, I^, we integrate 
the spin dependent differential conductance. The total 
current / = /| +/|, in units of /q = iet'^pt/nh, and the 
current polarization Pstm = cire shown in 

Figs. [2t and[2li, respectively. While the tunneling cur- 
rent / is small, the polarization Pstm can exceed 0.98. 
As we show below the possibility of injecting spin polar- 
ized currents is not restricted to the undoped case. 

In the case of doped graphene Ep is shifted from the 
Dirac point. In Fig. [3] we present results for the impurity 
spectral density A{u!) at low temperatures. The results 
show now a Kondo peak at Ep. Following Langreth [l^, 
it can be shown that the spectral density at Ep is given 
by the usual expression A{Ep) = sin'^ {nrid / 2) / nT where 
rid is the total charge displaced by the impurity and F = 
7Tp{Ep)V^y^. Our NRG results reproduce well this exact 
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FIG. 2: (Color online) (a) Impurity spectral density for the 
undoped case in the presence of a small magnetic field /iB = 
70jj,eV for the spin up (dashed line) and spin down (solid line) 
projections. Here Vhyb = lAeV, Ed ~ —U/2 and the other 
parameters are as in Fig. [l] (b) Same as (a) for Ed = — 0.5eV; 
(c) STM current and (d) current polarization as a function 
of the bias voltage Vt for the spectral densities shown in (a) 
(solid line) and (b) (dashed line). 



FIG. 3: (Color online) (a) Impurity spectral density for the 
doped case {Ef = — 0.35eV). The arrow indicates the position 
of the Dirac point, (b) Detail of the Kondo peak, (c) Impurity 
occupation as a function of the level energy; Viiyb = 1.4eV 
and Ef ^ -0.35eV (filled circles); Vhyb = 1.75eV and = 
0.35eV (filled triangles). The local interaction is in both cases 
U — 2eV. (d) Kondo temperature vs. Ed for the parameters 
of (c). The lines are parabolic fits. 



result. In Fig. [3j; we present results for the impurity 
charge Ud versus Sd- The impurity charge changes when 
the renormalized energies and 'sd + U cross Ep and the 
plateau corresponding to a localized spin in the impurity, 
Ud 1, is narrowed as Vhyb increases. This narrowing 
should not be interpreted as a reduction of the effective 
repulsion. As we show below the Kondo temperature is 
determined by the bare parameter U . 

We calculate the impurity spin susceptibility x{T), 
that at low temperatures shows universal behavior, and 
extract Tk using Wilson's criterion Tkx{Tk)/ IJ'^ — 
0.025. The results are shown in Fig. [SJi where log(Tx) 
as a function of Ed shows the usual quadratic behavior 
with a minimum at the center of the 1 plateau and 
a curvature determined by the bare interaction U . For a 
given set of impurity parameters, Tk varies with doping, 
or gate voltage, approaching zero for the undoped case. 

The transport properties of graphene are peculiar in 
many aspects and it is interesting to study the effect 
of adsorbed magnetic impurities. In what follows we 
present results for the resistivity pimp (T) due to these im- 
purities in the low impurity concentration (cimp) regime. 
Using the general expression for the conductivity in 
graphene [2C 



2l( and evaluating the band propagators 



Pimp{T) = PoVy 



hyb 



du J A{uj) 



duj 



(4) 



in the Born approximation we obtain, to first order in 



where po — nci^ph/e^. The general temperature de- 
pendence of the resistivity in the different regimes re- 
quires the numerical evaluation of the integral. As shown 
in Fig. m for the undoped case the resistivity is tem- 
perature independent while for the doped case we ob- 
tain the usual Kondo behavior. The NRG results can 
be qualitatively reproduced performing some simple ap- 
proximations. In the undoped case the spin depen- 
dent low frequency impurity spectral density is given by 
A^(w)~aV|^2yi, \u;\ [il-nda)/el+nds/ (ed+Uf] where Uda 
is the number of spin a electrons in the impurity orbital 
[23 |. For the sake of simplicity lets consider first the 
case Ed = —U/2 for which Uda = \, the resistivity is 
then given by pimp(r) = PaotV^y^ / e^. This result corre- 
sponds to impurities generating a short range potential 
of amplitude A oc 14iyb/k<i|- Nevertheless, in the ab- 
sence of electron-hole symmetry (e^ ^ —U/2) the system 
presents a large magnetoresistance due to the difference 
in the scattering rate of the two spin channels. The con- 
tribution of each spin is given by twice the r.h.s. of Eq. 
(HI) with A{lo) replaced by A„{bj)^ and the total resistivity 



Pi„,p (T, B) = pi,„p (T, B = 0) [1 - T^m^ (T, B)] , (5) 



4 




0.01 0.1 1 10 100 

T[°K] 



FIG. 4: (Color online) (a) Resistivity vs. temperature 
for doped {Ef ~ — 0.35eV, solid diamonds) and undoped 
graphene (solid squares) for B = 0, U = 2eV, Vhyh = 1.4eV, 
and Ed = — 0.63eV. Doped graphene shows Kondo scaling 
of the resistivity indicated by the dotted line, (b) Resis- 
tivity (triangles) and current polarization (circles) vs. tem- 
perature. Undoped graphene, jiB — 70^ieV (solid symbols) 
and = 700/ieV (open symbols). Other parameters are 
U = 3.5eV and ea = -0.5eV. 

with m{T,B) = {nd^-ndi) and for = 1, 7 = + 
U)^ -£d)/i{'=^d + U)'^ +8"^). The current polarization is 
P = {I^ — Ii)/ 1 = jm{T, B). These expressions are in 
good quaUtative agreement with the numerical results 
shown in Fig. 4b. The actual magnetoresistance and the 
degree of polarization of the current in bulk graphene 
will depend on the presence and intensity of additional 
scattering mechanisms. For the magnetoresistance to be 
observable, the scattering rate due to other mechanisms 
should be smaller than the one due to the impurities, that 
is, the impurities should be adsorbed on clean graphene. 

For the doped case the Kondo screening in the T ^ 
limit gives 

^ Cimp sin^(7rnrf/2) 
Pimp 0) = — , (6 

where n is the number of carriers per carbon atom. In 
the unitary limit, the resistivity is just determined by 
the ratio Ci^p/n. In the limit T ^ Tk we recover the 
resistivity characteristic of potential scattering defects. 
The temperature dependence of the resistivity, for T < 
Tk shows the universal Kondo behavior [23j . 

In summary we have numerically solved the problem 
of a magnetic impurity in graphene and analyzed the ef- 
fect of external in-plane magnetic fields. We find that 
as Ed varies, the region of stability for a localized spin 
{ud ~ 1) is narrowed and shifted with respect to the 
^lyb — * limit, in qualitative agreement with mean field 
results Kondo screening of the impurity spin oc- 

curs at low temperatures for the doped case. The Kondo 



temperature Tk decreases exponentially with decreasing 
doping and for the undoped case, the impurity behaves 
as a free spin down to zero temperature. We find very 
little Fano distortions in the STM spectrum that conse- 
cuently gives a direct measurement of the impurity spec- 
tral densities. Our central results concern the effect of 
magnetic fields: for zero or low doping, low Tk, the con- 
dition T,Tk < ^mB is accessible at moderate values of the 
external field B. In this regime the impurity spin is polar- 
ized and a non-magnetic STM tip can be gated to inject 
a spin polarized current, that is, the impurity acts as a 
spin valve. The magnetic field controls the degree and 
the axis of the spin polarization. In this regime, a finite 
impurity concentration leads to large magnetotransport 
effects in bulk graphene: for small values of \ed — Ep\ the 
system shows large magnetoresistance and spin polarized 
currents. All these effects are unique to graphene: they 
require l^y^^/|£d| to be large and the spin to be free at 
low temperatures, conditions that are never reached si- 
multaneoulsy in ordinary metals. Our main results are 
robust even in the presence of defects or other impurities 
that may change the structure of the pseudogap and the 
nature of the magnetic screening. However, for magnetic 
impurities close to one of these defects, new features are 
expected to appear in the STM spectra that will depend 
on the nature of the defect. 
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